SINGULAR COMPONENTS OF SPECTRAL MEASURES FOR 
ERGODIC JACOBI MATRICES 



C. A. MARX 

Abstract. For ergodic Id Jacobi operators we prove that the random singular 
components of any spectral measure are a.s. mutually disjoint as long as one 
restricts to the set of positive Lyapunov exponent. In the context of extended 
Harper's equation this yields the first rigorous proof of the Thouless' formula 
for the Lyapunov exponent in the dual regions. 



1. Introduction 

Given two sequences 

(1.1) (a n ) G /°°(Z; C) ,\a n \ > G Z , 

(1.2) (b n ) G r(Z;M) , 
we define 

(1.3) {Hip) n := h n i\) n + a n ^ n+l + a~ T^ n _i , 

a bounded self adjoint operator on T-L := / 2 (Z; C). 

Operators of the form ( 11.31) are known as 1 d - Jacobi operators (or Jacobi 
matrices), the most prominent and well studied case being Schrodinger operators 
where (a n ) = (1). 

We note that since f ll .3j) is unitarily equivalent to a Jacobi operator with a n > 
0, Vn, one could assume all a n > right from the start. In some situations it may 
however be more convenient to retain the original complex form, which is why we 
do not apply this unitary here. 

In this article, we shall be interested in an ergodic family of Jacobi operators, 
(H u ) uG q, induced by random variables a, b on a probability space (fl, J 7 , P), 

(1.4) a-.Q^C, b-.Q^R, 

(1.5) | IE^ log | a. 1 1 < oo , 

and an ergodic invertible map T on Q generating the random sequences (a n (uj)) 
and (b n (u))) according to 

(1.6) a n (u) := a(T n u) , b n (u) = b(T n u) . 

Remark 1.1. In particular, ( II. (jp together with (jl.5p imply that a n (u) ^ 0, Vn G Z, 
P-a.s. 
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We emphasize that in (II. 3D we do not require a positive lower bound for the 
sequence (a„). From a dynamical point of view, fixing z G C, this manifests itself 
in almost singularity of the 1-step transfer matrices 



Here, the phrase "almost singular" refers to the matrix a n B n whose determinant 
may come arbitrarily close to zero. For the ergodic counterpart in (jl.5p we even 
allow for the matrix a n B n (u) to be singular, however since a G L 1 this only 
happens on a set of zero probability. A more detailed discussion of singular Jacobi 
cocyles may be found in [TT] . 

As the main result of this article, we prove the following characterization of the 
singular(pp+sc)- spectrum for ergodic Jacobi matrices: 

Theorem 1.2. Let be an ergodic family of Jacobi operators satisfying U.5\) 
and U.6]) , and let d/x WiS be the singular component of some spectral measure for 
H w . Then, forF-a.e. uj and uj' , d/i W)S and d^Lj s are mutually singular on the set 
of positive Lyapunov exponent, {E : L(E) > 0}. 

Theorem 11.21 generalizes a result in |6j stated for the Schrodinger case. We note 
that this theorem will be proven using Kotani theory [TJ H], which usually gives 
results relating to the absolutely continuous component of the spectrum. It is also 
noteworthy that mutual singularity of the singular component of spectral measures 
is reminiscent of an analogous property exhibited by rank 1 perturbations of 
bounded s.a. operators (see e.g. [16]). 

We organize the paper as follows. Sec. [2] is devoted to the proof of Theorem ll.2l 
As a crucial first step, we prove that energies with positive Lyapunov exponent 
(LE) allow for two solutions of H^ip = Eip which are P-a.s. non-zero and 

I 2 (and hence necessarily exponentially decaying) at respectively +oo and — oo 
(Theorem 12 .4p . 

Consequently we obtain that given an energy E with positive LE, the probabil- 
ity for E to be in "the natural" supporting set (defined in (I2.19P ) of the singular 
component of any spectral measure is zero (Theorem 12.61) . This result is known 
for Schrodinger operators [I]. 

While following the basic strategy laid out in [I], our proof even when applied 
to the Schrodinger case differs in a number of technical details. In particular, 
Theorem 12 .41 strengthens the respective statement for Schrodinger operators given 
in Theorem 6.5 of [I], by removing a set of zero probability, thus simplifying 
subsequent arguments (for details see remark |2"3|) . 

Finally, Theorem 11.21 follows as a corollary of Theorem 12.61 

Sec. [3] is devoted to an application of Theorem 11.21 which really served as 
motivation for the results of Sec. [21 Ergodic Jacobi operators with singular 
transfer matrices arise naturally in a tight binding description of solids subject to 
an external magnetic field. 



(1.7) 
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A prominent example is extended Harper's model [121 El 13 [TH] , associated with 
the following almost-periodic operator Hj e . a defined on Z 2 (Z): 

(1.8) {H Xe .J)) k := v(9 + ak)i) k + c(0 + ak)ip k+1 + c(6 + a(k - l))^ fc -i , 
where A = (Ai, A2, A3) and 

(1.9) cj(x) := \ 3 e~ 2n < x+ ^ + A 2 + Aie^+f) , v(x) := 2cos(27rs) , 

are analytic functions on R/Z. 

The operator in (11.81) describes the influence of a transversal magnetic field to 
the motion of an electron in a 2-dimensional crystal layer. Here, a is a fixed 
irrational number associated with the magnetic flux through the unit cell and 
9 G [0, 1) is the (random) quasi- momentum usually referred to as "random phase". 

The parameter triple (Ai, A2, A3) {coupling) models the lattice geometry as well 
as interactions between the lattice sites allowing for both nearest (A2) and next- 
nearest neighbor (Ai and A3) interactions. Without loss of generality one may 
assume < A2 , < A1+A3 and at least one of Ai, A2, A3 to be positive [9]. We note 
that the quasi-periodic operator in (I1.8P specializes to the famous almost Mathieu 
operator (in physics literature also known as Harper's model) upon setting Ai = 
A 3 = 0. 

Even though some attempts, both rigorous [121 E] an d heuristic [3 HE] , have 
been made, so far little is known about extended Harper's model. As an important 
ingredient for the spectral analysis it is desirable to obtain the Lyapunov exponent 
(on the spectrum) as a function of the coupling constants. 

Theorem 11.21 enables us to prove a formula of the Lyapunov exponent in the so 
called "dual regions" of parameter space (Theorem 13. ip . A precise meaning of the 
word "duality" is given at the beginning of Sec. [3j We mention that Theorem 13.11 
provides a rigorous proof of results obtained in physics literature by Thouless [7J . 

Finally, we mention that the missing link to a complete understanding of the 
Lyapunov exponent, the so called "self dual" region of coupling constants, a priori 
cannot be approached via duality. Computing the Lyapunov exponent in the 
self dual region therefore constitutes an important open problem which so far 
has escaped treatment even on a heuristic level in physics literature [7J. The 
development of these a priori non-duality based methods will be the subject of 
future work. 

Acknowledgement. I would like to thank Svetlana Jitomirskaya for valu- 
able discussions during the preparation of this manuscript. I would also like to 
thank the anonymous referee for insightful remarks leading to improvement of our 
presentation. 

2. Proof of the Main Theorem 

We start with the following basic result valid for any ergodic family of bounded 
self adjoint operators: 
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Theorem 2.1. For fixed d G N, let H w be any ergodic family of bounded self 
adjoint operators on Z 2 (Z d ). Assume that for F-a.e. u, the eigenvalues of H u 
have finite multiplicity, whenever H w has non-empty point spectrum. Then, fixing 
E el, F{u : E is eigenvalue for H^} = 0. 

Proof. Let B G R be any fixed Borel set and denote by Pb(H w ) the spectral pro- 
jection onto B. The claim follows since by ergodicity dimRanP^i^) is constant 
P-a.s with value either oo or [HI [3]. □ 

In order to relate the spectral properties of the family H u to its dynamics given 
in terms of the Lyapunov exponent, we shall consider solutions u(z) in C z of the 
second order difference equation Hu = zu, where z G C is fixed. For two such 
solutions u and v (same z), define the Wronskian 

(2.1) Wntu,v] := u n v n+1 - u n+ iv n ,n G Z . 

It is easy to see that the Wronskian of two solutions satisfies the following con- 
servation law 

(2.2) a^~[W n -i[u, v] = a n W n [u,v] ,n G Z . 
As an immediate consequence of (12.21) we obtain 

Proposition 2.1. Any eigenvalue of a Jacobi operator has multiplicity one. 

We note that Proposition 12.11 uses a n ^ 0, Vn G Z. In particular, in the ergodic 
setup, remark [1J] implies that Proposition 12. II holds P-a.s. Simplicity of the point 
spectrum of Jacobi operators shall be crucial for Sec. El 

Sometimes it is useful to relate the solutions of the Schrodinger equation for 
the double-sided and the respective one-sided operator H + := P + HP + , where P + 
is a projection onto the positive half line. 

For z G C fixed, a solution if) to Hip = zip is a solution of the respective 
(positive) half-line equation if and only if it satisfies the boundary condition ip = 
0. In particular this implies that Proposition 12.11 also applies for eigenvalues of 
H+. 

More generally, for any n G Z one can consider the positive n-one-sided operator 
H™ := p+ n HP +n with P +n the projection onto {m > n}. For convenience, we 
allow n = — oo in which case we identify H™ with the double-sided operator 
H. Similarly, one defines the negative n-one-sided operator, := P~ n HP~ n , 
n G Z, where P~ n is the projection onto {m < n}. Again, we adopt the convention 
H = H™ for n = +oo. To simplify notation we write H± := H±. 

For the ergodic setup, Proposition 12.11 and Theorem 12.11 imply the following 
statement which will be key for the proof of Theorem 11.21 

Proposition 2.2. Let H u be a family of ergodic Jacobi operators, and let E G K 
be fixed. Then the probability for E to be an eigenvalue of some (positive or 
negative) n-one-sided operator (including n = ±oo) is zero. 
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For z G EI + , let if>+(z) be a solution to the equation Hip = zip which is I 2 
at +00. Self adjointness of the half-line operator implies existence of ip+(z) as 
well as ip +>n (z) 7^ 0, Vn. Using (12.21) . we conclude uniqueness of ip+(z) up to a 
multiplicative constant. Similarly, we can find a solution ip-(z) which is I 2 at —00. 

It is well known that the Green's function of H can be expressed in terms of 
ip ± (z) by 

2.3 G(n, m; z) := (8 n , G{z)5 m ) = , 

a m W m [ip^{z),ip + (z)\ 

for z e C \ K. Here, {5 n , n£Z} denotes the standard basis in Z 2 (Z). 
Finally, for n 6 Z we define the m-functions, 

(2.4) m + ,„ W := -%±44 , m., n( z) := -J^lW , 

where we set m_|-(2;) := m± t Q(z). It is clear that in the ergodic case we have 
m ±,n(w) = m±(T n u), n G Z. For simplicity of notation, we will suppress the 
z-dependence of the m-functions whenever the context permits. 

Since ip± are solutions to the Schrodinger equation, the m-functions satisfy the 
following Riccati-type equations 

h (b n - z) - \a n \ 2 m +<n = 

(2.5) -|an-i| 2 ^-,n + (&n - z) — = . 

"7.-,n+l 

We mention that using a relation analogous to (12. 3p (see e.g. [3 E]), the 
functions m±(z) relate to Borel transforms of the spectral measures dfi± associated 
respectively with H± and 8±i, i.e. 

(2.6) m±(2) 



.r 



This in particular implies that m± are Herglotz functions, and we may use the 
theory of Borel transforms of measures to consider limits as Im(^) — > + . 

From here on we specifically consider the ergodic set-up. The following state- 
ment relates the Lyapunov exponent and the m-functions. It extends a result 
known for the Schrodinger case |15j . 



Proposition 2.3. Forlmz > we have 
(2.7) 2L(z) =E w log(l+ Jm: 



|a (u;)| 2 Imm + (w, z) 



Remark 2.2. We mention that a formula similar to (I2.12p is given in [T7] (see Eq. 
(5.73) therein), however with |a (o;)| 2 in (12.31) being replaced by a (u). Since a 
proof is omitted in [T7] , we correct the typo in [T7] providing the details below. 
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Proof. Taking imaginary parts, the first equation in (12. 5ft implies 



log <^ \a n (u)\ 2 + 



Im z 



logImm + (T™ u 



Im m + (T n uj) 
-loglm m + (T n u) - 2 log|m + (T n - 1 a;)| 

Since T is measure preserving, taking expectations on both sides yields 



(2- 



e„ log a + 



Imz 



a (a;)| 2 Im m + (uj 
First, notice that by Birkhoff's ergodic theorem 



2E W log 



(2.9) 



71 



log 



-> ^ lOg 



VM-,o(w) 



P-a.e. , 



as n — )■ oo. 

On the other hand Imz > implies L(z) > 0, hence using Oseledets' theorem 
(in its original GL(2,C) formulation [13]) we obtain 



(2.10) 



n 



log 



-L , P-a.e. 



as n — > oo. Thus, equations (12. 8p . (12. 9p . and (I2.10p finally yield the claim. □ 
Using ( 12. 6p . the monotone convergence theorem implies the existence of 

(2.11) S±(u,,E):= lim + , 

for any fixed Bel. 

Thus using dominated convergence, the formula for the Lyapunov exponent 
given in Proposition 12.31 persists in the limit, i.e. 



(2.12) 



2L(E) =E u ]og<l + 



\a (uj)\ 



; S + (u,E) 



-i 



for any fixed E G R. Here, we also used that for any E G R, LiE + ie) is 
continuous w.r.t. e which e.g. can be obtained from the Thouless formula. 



Since (12.121) contains the limit of the quantity 
will prove to be of use: 

Lemma 2.3. For Imz > ; 

Imm + (w,z) 1 

(2.13) 



Im z 



the following Lemma 



Imz 



\a (uj)\ 2 \ip + ,o(uJ,z)\< 



\P+if> + (u,z) 



Proof. Using functional calculus, Eq. (12.61) implies 
Imm+(w, z) 

n=l 



(2.14) 



Imz 



J]|<5 1 ,(P + i/ w P + -z)- 1 5 n )| 
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Finally, making use of expressions for the Green's function of the half line operator 
(see e.g. [15J), we obtain 

for any n > 1. □ 



(2.15) \(6 1 ,(P+H u P+-z)- 1 5 n )\ 



We are now in the position to characterize the set of energies with positive 
Lyapunov exponent in terms of the m-functions. 

Theorem 2.4. Let E G M. s.t. L(E) > 0. Then for P-a.e. uj, there exists a 
solution u+(u),E) with u+ jn (u,E) ^ for all n G Z, which is I 2 at +oo. In 
particular, m + (u,E + iO) exists, is finite and non-zero. 

Remark 2.5. In [3], the authors proof a similar statement for Schrodinger opera- 
tors, which however does not rule out the case where the limiting solution u + (u, E) 
has zeros (Theorem 6.5 (b) [3]). We mention that even though the authors of [3] 
suspected this to only happen with zero probability, no proof was provided then. 

Proof. Using (12.121) . L(E) > implies that S + (oj,E) < oo on a set of positive 
measure P; in particular, letting Q := {u : S + (T n u, E) < oo for some n G Z} we 
obtain P(fi Q ) > 0. In fact, since T _1 f2 = Q , we have P(fio) = 1 using ergodicity 
of T. 

Let U) G Qq and n£Z such that S+(T n u, E) < oo. We mention that for n = 0, 
we already obtain the statement using a standard argument from the the theory 
of Borel transforms. 

Using Lemma [2.31 for z = E + it and e > 0, 



Imm + (T"u,z) 1 

2-16 f ^ = r, PA (T n u,z 2 , 

Ini2 |ao(j n w)| 



where we define V ; + to satisfy the condition ip +!0 (T n u, z) := 1. 

Then, by Lemma 12.31 and Banach Alaoglu (S + (T n u, E) < oo), we obtain a 
limiting solution u+(u,E) of the Schrodinger equation, 

(2.17) w - lim ip + {T n u, E + ie m ) =: u+(u, E) , 

m— >oo 

for some sequence (e m ), e m \ 0. By construction u+(u,E) is / 2 at +oo. Making 
use of f)2.2p . any solution of the Schrodinger equation which is I 2 at +oo is uniquely 
determined up to a multiplicative constant. Since by construction u + {n) = 1, this 
limiting solution is in fact independent of the sequence (e m ), whence 

(2.18) u+(u, E) = w- lim ip + (T n u, E + it) . 

If u+ jm (oo, E) = for some m G Z, -E) is an eigenfunction of H™ u . By 

Proposition 12.21 this however happens with zero probability. 

Thus, P-a.s. we have u +t7n (u, E) ^ 0, Vm G Z. In particular, taking the limit 
in the definition of m+, we obtain the claim. □ 
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Following Deift and Simon [3], for u G we consider the following set 
(2.19) 

6 W := \e g R : lim sup {Im G w (0, 0; £ + ie) + ImG w (l, 1; E + ie)} = oo 

It is well known that the singular component of any spectral measure for H u is 
supported on & w (Theorem of de la Vallee-Poussin; see also [17]). 
G(0, 0) and G(l, 1) relate to the m-functions by 

— G(0, 0; z)~ l = \a \ 2 m + (z) + \a^i\ 2 m_(z) + z — b 

(2.20) 67(1, l;*) -1 = m+iz)- 1 + \a \ 2 {\a„ 1 \ 2 m_(z) + z-b }- 1 . 

These relations are obtained directly from the expression for the Green's function 
given in (12.31) . 

The following statement extends the statement of Proposition 12.21 to <S W , the 
natural supporting set of the singular component of any spectral measure for 
H u . The respective result was first proven for Schrodinger operators by Deift and 
Simon [1]. 

Theorem 2.6. Let E Gi such that L(E) > 0. Then, ¥{oo : E G ©4 = 0. 

Proof. By Theorem 12 A\ we know that for P-a.e. u, the limits m±(uj,E + i0) both 
exist, are finite and non-zero. Using the relations (12.201) for such u, E 6 & u if 
and only if at least one of the following happens 

\a (u)\ 2 m + (u,E + i0) + \a_ l (u)\ 2 m_(uj,E + i0) + E -b (u) = 0, 

m + (u,E + + \a (u)\ 2 {\a^ 1 {u)\ 2 m_{u,E + i0) + E -boiu)}- 1 = 0. 

Using the definition of m±, it can easily be verified that both these situations 
imply that the respective limiting solutions u±(u,E) constructed in the proof 
of Theorem 12.41 are in fact constant multiples of one another. Thus we obtain a 
solution which is I 2 at both ±oo, i.e. an eigenvector of H u . Employing Proposition 
I2.2[ the set of such u has zero probability, which yields the claim. □ 

Finally we obtain Theorem 11.21 which concludes this section. 

Proof of Theorem ! Using the definition of the Lyapunov exponent, the set 
£:={£Gl:0< L(E) < oo} is measurable. Consider, W. := {(u, E) e Q x £ : E G 
Letting f(u,z) := ImG^(0,0;;z) +lmG u (l,l;z), where (u,z) eQx (£ + iR+), 
we note that 9T can be written as 9t = {(co, E) G VL x £ : lim sup e _ i>0 + f(co, E+ie) = 
oo}. Since / is measurable, we conclude measurability of 01. 
The rest now follows using Tonelli. By Theorem 12.61 Vu/ G Q : 

= [f{lo:Ee &M^> 

= y^(£ne w )dP(a;) , 

which implies the claim. □ 
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3. AN APPLICATION: DUALITY AND EXTENDED HARPER'S MODEL 

In this section we illustrate the use of Theorem 11.21 to obtain an expression for 
the Lyapunov exponent in the dual regime of extended Harper's model . The 
model as well as its relevance has been discussed in Sec. [TJ 

Similar to the almost Mathieu operator, the Hamiltonian in (jl.8p exhibits an 
intrinsic symmetry in the coupling constants under Fourier transform, known as 
duality. Heuristically speaking, duality relates certain regions in parameter space, 
mapping localized states to Bloch waves. 

To give a precise meaning of duality, on the Hilbert space W := L 2 ([0, 1) x Z), 
consider the unitary operator U defined by 

(3.1) {U(f)){r],m) := X) / d6e i2nmd e i27Tn{ma+v) 0(0, n) . 

ngZ 

We note that such operator has originally been introduced by Chulaevsky-Delyon 
[2] and later employed in [6], both in context of the almost Mathieu operator. 

For A 2 7^ 0, let a denote the map, <r(A) := ^(A 3 ,l,Ai), where as earlier A = 
(Ai, A 2 , A3) represents the coupling in extended Harper's model (see (jl.9p ). 

We also define the following regions in parameter space: 

region I: < Ai + A 3 < 1, < A 2 < 1 , 
region II: < Ai + A 3 < A 2 , 1 < A 2 , 
region III: max{l, A 2 } < Ai + A3 . 

This partitioning is induced by the map a since 

(3.2) o[l°) = 11° , a(III) = III . 

In this sense, the interiors of region I and II are dual, whereas region III is self-dual. 
For a more general discussion of duality we refer to |12j . 

In [9] (see Theorem 1 therein) it was proven that for Diophantine a and A G 1°, 
the spectrum of Hj g . a is only pure point (i.e. for a certain set T of phases 9 
with I To I = 1). Recall, a G 1 is called Diophantine if there exists < b(a) and 
1 < r(a) < +00 s.t. for all j e Z \ {0} 

(3.3) |sin(27rj«)| > b ^ 



\j\r(a) 



It is known that DC(r) := {a : (13. 3p holds for r and some b} is a set of full 
Lebesgue measure for every r > 1. 

As we will argue, Theorem 11.21 will allow to relate the regions I and II estab- 
lishing the following: 



Theorem 3.1. For every Diophantine a, the LE in region II is zero on the spec- 
trum. In region I, the Lyapunov exponent on the spectrum is positive and given 
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by the formula, 



1 + Vl-4AiA 3 




2Ai 



if Ai > A 3 , A 2 < A 3 + Ai , 



(3.4) L{a; A, E) = < log ( 1 + — ^- 4AlAa ^ , 1/ A 3 > A x , A 2 < A 3 + Ai , 



1 + a/1 - 4AiA 2 



A 2 + v^l - 4A i A s 



«/ A 2 > A 3 + Ai . 



H in (11. 8p naturally induces a bounded self adjoint operator on H.', which shall 
be denoted by Hi , 

J A;cr 

(3.5) HLJ:=H- X . 6i j(9,.) , V0 G [0, 1) . 

Using the definition of the unitary [/, direct computation establishes the fol- 
lowing relation between dual regions: 

(3-6) UH'^U- 1 = A 2 #; (X);Q • 

To avoid confusion, we emphasize that (13.61) shows that Hi and H' ^ are 

unitarily equivalent and not the respective operators on the smaller space H = 
l\Z). 

An important property of duality is that preserves the density of states [i~2] . 
To see this in the present formulation notice that 

(3.7) U8 = 5 , 

where 5q is viewed as element of H', i.e. employing the embedding H^H! 

(3.8) (u n ) 1 — y {<p{6, n)) , where <f>(6, n) = u n e Z , 9 G [0, 1) . 

Denote by dn({Hj. e Q } ; E) the density of states of the ergodic family {Hj. e Q | 
defined as usual by 

(3.9) Eg <5 , f(H- XAa )5 ) =: J f(E)dn({Hj. e J ; E) , 

for / continuous and compactly supported on M. Then, combining (I3.6p and ( 13. 7p . 
we obtain the desired invariance of the density of states, 

(3.10) dn {{H Xe J ;E)=dn ({ \2H a ^ a ) ; e) = dn ({fl^^} ; A,- 1 ^) . 

We mention that the last equality in (13.101) follows from the spectral mapping 
theorem. 

In [6] it is shown that fixing A and a, it is possible to find a labeling of the 
eigenfunctions of H xe . a that produces measurable functions w.r.t. the random 
phase 9. 

The proof in [6] (Theorem 2.2) is given for an arbitrary bounded, ergodic 
Schrodinger operator on l 2 (1i), the crucial property needed being however only 
simplicity of the point spectrum. As shown in Proposition 12.11 the same is true 
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for Jacobi operators. Therefore, Theorem 2.2. in [6] extends to an ergodic family 
H u of Jacobi operators, as defined in Sec. [2J 

A measurable labeling of eigenfunctions is obtained as follows: For u G Q, let 
u be a normalized eigenvector of H u . Pick the unique j G Z such that \u(j)\ > 
\u(k)\, VA; and \u(j)\ > \u(k)\, V/c < j ("left-most maximum"). We say u is 
attached to j. Let Nj(u>) be the number of eigenfunctions attached to j, and 
for j, k G Z set Vtj^ '■= '■ > k}. On f^i, let (ui(n; u, j)) be the 

eigenfunction attached to j with maximal jth component (in magnitude) and 
ei(u,j) its associated eigenvalue. If multiple eigenfunctions attached to j have 
the same magnitude of their jth entry, discriminate them by choosing the one 
with largest eigenvalue. Note that since the point spectrum is simple, there can 
only be finitely many different eigenfunctions attached to j whose jth components 
coincide in magnitude. For further details we refer to [6]. 

For this labeling of eigenfunction we have (see Sec. 2 in [6]), 

Theorem 3.2. Let be an ergodic family of Jacobi operators defined as in Sec. 
[H and Nj(u)),ei(oj,j), and (ui(n;u,j)) n( zz defined as above. Then, fixing l,j and 
n, Nj(u),ei(u,j), and ui(n; u) , j) are measurable functions on Q. 

The next statement shows that applying the unitary U defined in (I3.ip . we 
obtain spectral measures that are constant for a.e. phase. The proof carries over 
from [5] (see Theorem 3.4 therein) without any modification. 

Theorem 3.3. Fix A and a Diophantine and let (ui(n;9,j)) be the measurable 
labeling of eigenfunctions for H a( j^ e . a described in Theorem \3.2L For f{6) G 
L 2 ([0, 1], d9) arbitrary and some fixed I, j G Z, set (f)(0,n) := f(9)ui(n; 9, j) . De- 
note by dfirj the spectral measure for Hj g . a and ip^n) := (U(j))(r],n) . Then, dfi^ 
is i] independent for a.e. r\. 

Equipped with Theorem 13.31 we are now ready to prove Theorem 13.11 

Proof of Theorem \3.1\ For a fixed Diophantine frequency a, let £± ;= {E G spec : 

L(a; A, E) > 0}. We claim that £± = for any A in region II. 

First consider A G IP and let d^ be defined according to Theorem 13.31 Then, 
since dfi v is constant a.e. rj, Theorem 11.21 implies that in fact d^, v is purely a.e. 
on £± for a.e. rj. Hence, using a well known argument by Pastur-Ishii (JSJ [T4] ; 

see also the proof of Theorem 9.13 in [3]) we obtain P^ c (£i) = P^(£±) = 0; here 
P V (B) and P^ C (B) denotes the spectral projection associated with a Borel set B 
and its restriction to the ac subspace of H v , respectively. In particular, in terms 
of the density of states this implies 

spec = supp {dn {{H J e . a }) } C I \ £+ , 

which by continuity of the Lyapunov exponent w.r.t. to E [HI [10] yields the claim 
for A G IP. To obtain the claim for all off region II, i.e. including the boundary, 
we employ the following continuity argument: 
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Let A be an arbitrary point in dll and E a fixed energy in the spectrum of 
Hj. Using continuity of the spectrum in the Hausdorff metric 0, there exists a 
sequence (A n ), approximating A from within 11°, and a corresponding sequence 
of energies (E n ), with E n in the spectrum of Hj , such that E n — > E. Since the 
Lyapunov exponent is jointly continuous at (E, A) for a fixed Diophatine a [TTj . 
we obtain L(a; A, E) = 0, as claimed. 

Making use of Thouless' formula and (13.101) , zero Lyapunov exponent in region 
II immediately implies an expression for the dual region since 



L(a; A, E) = - J log\c x (x)\dx + J log\E - E'\dn({Hj Aa };E') 



(3.11) 



log 



(A) I 



dx + L(a;a(X),\2 1 E) 



Here, cj is defined as in (11.91) . 

Using Jensen's formula, the integral in (13.1 ip can be computed explicitly [9] 
(see Sec. 2, p. 107), 
(3.12) 



log|cr(a;)|dx 



'log A 3 
logAi 

log 



if A 3 > Ai > 

and Ai + A 3 > A 2 > , 
if Ai > A 3 > 

and Ai + A 3 > A 2 > , 



2A!A 3 



-A 2 + y/Xl - 4AiA 3 
giving rise to the expression for the Lyapunov exponent in region I given in ([3 



if Ai + A 3 < A 2 and Ai, A 3 ^ 
if Ai + A 3 < A 2 , Ai or A 3 = . 



□ 
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